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We find an algorithm of numerical renormalization group for spin chain models. The essence of this algorithm 
is orthogonal transformation of basis states, which is useful for reducing the number of relevant basis states to 
create effective Hamiltonian. We define two types of rotations and combine them to create appropriate orthogonal 
transformation. 



1. INTRODUCTION 

A method of Hamiltonian diagonalization is 
suitable for calculating physical quantities associ- 
ated with wavefunctions, such as structure func- 
tions and form factors. However, we cannot diag- 
onalizc Hamiltonian directly in quantum field the- 
ories because dimension of Hilbert space is gener- 
ally infinite. To create effective Hamiltonian, we 
extend a technique of NRG (numerical renormal- 
ization group) proposed by K. Wilson p^. 

In spin chain models, any state can be ex- 
panded with a direct product of two sets of ba- 
sis states, each of which is for one of two spin 
blocks. For a finite lattice, we can calculate wave- 
function of the ground state by diagonalizing 
Hamiltonian. In Ref. j^j, S.White has found that 
singular values D k of target-state wavefunction 
i&ij = J2k UikDkVkj are useful for the purpose of 
creating effective Hamiltonian. 

i j k 

where \u k ) = Y,i U ik\i) and \v k ) = Y,j V kj\j)- If 
|-Dfe| is large, the corresponding basis state is rel- 
evant. However, if \D k \ is small, the correspond- 
ing basis state is not relevant and can be thrown 
away from calculation. We can control calcula- 
tion accuracy using singular values D k when we 
truncate Hilbert space. Magnitude of \D k \ deter- 
mines the relevance of basis states in the target 
state. Basis truncation based on singular values 
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works well in one-dimensional spin models such 
as Heisenberg and Hubbard models |2I3| . This is 
an extension of Wilson's NRG and called DMRG 
(density matrix renormalization group). White 
has shown that there exists orthogonal transfor- 
mation that largely reduces the number of basis 
states with high calculation accuracy. However, 
we cannot assure that a number of singular val- 
ues take values of nearly zero in general cases. In 
order to create optimized basis states without re- 
lying on density matrix, we will find another use- 
ful basis transformation. Our aim is to find basis 
transformation that make coefficients (i.e. trans- 
formed wavefunctions) exactly zero when target 
states are expanded with new basis states. 



2. FORMULATION 

Let us consider S = 1/2 Heisenberg chain. 
Since we cannot diagonalize infinite dimensional 
Hamiltonian directly, we start from a finite lat- 
tice. We extend a lattice step by step by adding 
one site to it in each RG transformation. \ui) 
is a complete set of basis states for a finite lat- 
tice L. Hilbert space of an extended lattice with 
size (L + 1) is spanned with a set of basis states 
|u<)|j}j where |j) is an eigenstate of z-component 
of a spin operator for the (L + l)-th site 

S z \j)=±±\j) for j = ±. 
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We want to find a rotation matrix U that decrease 
the number of relevant basis states. 



(2) 



In general, any orthogonal matrix can be ex- 
pressed as a product of two-dimensional rotation 
matrices R. 
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This matrix rotates the fc-th and Z-th basis states 
by an angle 9. Let us consider a ground state that 
is an eigenstate of a Hamiltonian. 

|*> = 5>«l«i>|j>- (3) 

ij 

We consider a part of this target state. 

|*}(i,2) = *i-|ui>|-> + *i+|ui>|+> 
+*2-|«2>|-) + *2+|u 2 >|+>. 
If we rotate two basis states and \u 2 ) 

\u\) — cos9\ui) — sin 9\u' 2 ), 
\u 2 ) =sm9\u[) +cos9\u' 2 ), 

we have 

l*)(l,2) = *i-|«i>|-> + *i + |«i>|+> 
+ *2-l«2>|-> + *2+l«2>l+>, 

where W are transformed wavefunctions 

= cos 9 + * 2 - sin 6, 

= *i+ cos 9 + * 2 + sin 9, 
* 2 _ = sin 9 + <3> 2 - cos 6, 

* 2+ = sin 9 + <3/ 2 + cos 9. 

We define the following two types of rotations. 
Rotation I is defined to make the fourth compo- 
nent zero 4 , 2 + = 0; 

*i+ I i | *i + 
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A condition ^' 2 + = determines the value of 9. 
Rotation II is defined to make the third compo- 
nent zero ^2- = when components and 
^2+ are zero; 
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When Wi+ = $> 2 + = 0, the transformed compo- 
nents and \D r 2 + are always zero. Therefore, 
we can make ^> 2 _ zero by choosing 9 appropri- 
ately. 

In order to understand how these two rotations 
are used to find relevant basis states, let us con- 
sider a simple example where a target state is 
expanded with six basis states. 
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The following is transformation of the six- 
dimensional vector of the target-state wavefunc- 
tion i&ij. 
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We obtain the final form by transforming the 
vector with two rotations appropriately. If a 
transformed wavefunction is zero, the correspond- 
ing basis state is irrelevant and can be ne- 
glected. That is, only the first three basis 
states (|m")|±), |m 2 ')I — )) are relevant and oth- 
ers (|u 2 ')|+), | it?) } | ±)) are redundant. The for- 
mer three relevant basis states can reproduce the 
target state exactly without other basis states. 
There exists orthogonal transformation that re- 
duces the number of relevant basis states into 
three without causing calculation error. This is 
also true for arbitrary dimensions other than six 
of this example. 

Let us consider a case where we target ground 
and first-excited states. The followings are wave- 
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function vectors of the target states. 
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The indices (1) and (2) indicate the ground and 
first-excited states, respectively. We optimize ba- 
sis states in two steps. In the first step, we target 
the ground state \E'( 1 ) and do not take care of 
the first excited state f™. In the same way as 
the previous example we can find three rele- 
vant basis states that can reproduce the ground 
state exactly by using a combination of rotations 
I and II. In the second step, we target the first ex- 
cited state to optimize basis states further. 
When we rotate basis states, we don't touch the 
first four basis states (|i4)|±), |it2)|±)). We ap- 
ply a combination of rotations I and II to other 
(2N — 4) pieces of basis states ((1*3) |±), It/4) j±), 
. . ., |u^)|±)). Finally, we obtain seven basis 
states (K)|±), K)l±>, K')l±>, lOI-)) that 
are relevant for reproducing the ground and first 
excited states exactly. In the same way, we can 
also create relevant basis states when arbitrary 
number of eigenstates are targeted. The number 
of target states can be chosen according to the de- 
sired calculation accuracy. We can use this tech- 
nique of basis transformation to truncate Hubert 
space and create effective Hamiltonian in NRG. 
Our NRG algorithm for S = 1/2 Heisenberg chain 
is as follows: 

1. Diagonalize Hl {Hlu u = E n u n ), which is 
Hamiltonian for a finite lattice L. 

2. Transform basis states u n — > u' n using type 
I and II rotations to make wavefunctions of 
target states zero as many as possible. 



3. Create effective Hamiltonian 
0{H L L , where L = 
Dimension of Hl is M. 
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4. Add one site to Hl to provide next Hamil- 
tonian Hl+i- Dimension of Hl+i is 2M. 

5. Regard Hl+i as Hl, and go to step 1. 

These steps are iterated till vacuum energy per 
site converges with the desired accuracy. I have 
applied this algorithm to S = 1/2 Heisenberg 
chain and calculated vacuum energy per site. An 
exact value of vacuum energy has been obtained 
using Bethe ansatz. The difference between our 
result and exact value is about 3.9 x for pa- 
rameters M = 256 and L = 2000. 

3. CONCLUSION 

We have found basis transformation that 
largely reduces the number of relevant basis 
states. Especially, three states are sufficient for 
reproducing one target state in each RG step. 
Relevant basis states reproduce target states ex- 
actly in each RG step. Therefore, if dimension 
of effective Hamiltonian M is sufficiently large, 
calculation errors come from non-targeted states. 
Future problems are improvement of accuracy, 
calculation of correlation functions, calculation of 
Haldane gap, and application to other models. 
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